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ABSTRACT 


In a space booster on takeoff, a control system must be 
employed to prevent the rocket from falling over as it is 
forced upward by the engines. One accurate dynamic model of 
the space booster on takeoff is the inverted pendulum. 

This paper investigated the inverted pendulum system from 
development of equations of motion, through implementation of 
an actual inverted pendulum system using state variable feed- 
back control. The concept of state variable feedback was 
analyzed in determining a solution to the problem of control 
of the inverted pendulum. 

Most importantly, design and construction of a working, 
solid state controlled, inverted pendulum demonstration model 


was accomplished. 
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IT. INDRODUCTION 


Mipememsevay Of Automatic Control Systems, it is not 
uncommon to encounter a plant which is inherently unstable. 
This paper investigated the inverted pendulum as an example 
of an unstable system which could demonstrate many of the basic 
aspects of control theory, and at the same time be constructed 
for use as a classroom model in Naval Postgraduate School's 
Aeronautical Engineering Advanced Control Theory course. 

The inverted pendulum is a model of the attitude control 
of a space booster on take-off. The objective of the atti- 
tude control is to keep the space booster in a vertical 
Becwtion. 

The system considered was a four-wheeled cart on which a 
bearing mounted inverted pendulum was installed. For this 
thesis, a cart was designed and constructed keeping in mind 
that overall size and weight should be minimized by using 
current solid state technology for the electronic controls and 
power supply. 

The cart is driven by a small D.C. servomotor which is 
geared to one pair of wheels through an axle. The cart can 
move in only one horizontal axis in both positive and negative 
directions. The voltage applied to the motor comes from sen- 
sor Signals describing the position and motion of the free 
pendulum and the position and motion of the cart. 

The sensors used are a permanent magnet tachometer geared 


to the pendulum to measure rate of pendulum movement about the 





pivot, potentiometers which measure pendulum angle and cart 
position (using the second axle/wheel combination), anda 
techometer within the motor housing to measure cart velocity. 

The pendulum/cart assembly represents the plant, and the 
plant is unstable since the pendulum cannot remain upright 
without help of a force from the motor/wheel combination to 
keep the cart directly underneath the vertical pendulum. The 
objective of the automatic control system for the inverted 
pendulum system is to keep the pendulum in as nearly a verti- 
cal position as possible while returning the cart itself to 
mes Starting position. 

Analysis of the system started with the development of the 
equations of motion for the inverted pendulum and cart; ex- 
pression of the linearized and simplified equations in state 
variable notation; development of a control law for full 
state variable feedback; and construction of a working model 


of the inverted pendulum. 





fio, PROeGEM FORMULATION 


fae UNVERTED PENDULUM EQUATIONS OF MOTION 

Consider an inverted pendulum with the pivot of the 
pendulum mounted on a motorized cart which can move in the 
MOrizontal direction. The cart is drive by a small D.C. 


Servomotor that exerts a force u(t) on the cart. 


Speeem Variable Notation (refer to figure 1) 


s(t) - displacement of the pivot point 

9(t) - angular rotation of the pendulum from vertical. 

m - mass of the pendulum 

M = mass of the cart 

2 - total length of the pendulum 

L - distance from pivot to center of gravity of the 
pendulum 

a - effective pendulum length. (L' = J+mL") 

J - moment of inertia of the pendulum with respect to 
phe c.g. 

H(t) - horizontal reaction force in the pivot 

V(t) - vertical reaction force in the pivot 

u(t) - force exerted on the cart by the motor 

F - viscous friction coefficient of the cart (includes 


friction, gearing losses, back EMF of the motor, etc.) 


g - gravitational acceleration 


Pendulum Equations 


2 
Tee ee ein 8(t)) = «8 (t) (1) 
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)F = m oe IL Ceeema ce); = V(t) = mg (2) 
Y dt 
a*s (t) 
eee it) Gsin g(t) -H(t) Leos 9(t) () 
ele 


eact Equations 


ds (t) 
dt 


Blogs eo) 


at? 


u(t) - H(t) - F 





= = OM 
= 5 
Be, 0 (5) 


iemertctren tem F 1s due to the cart motion. The fric- 
tions of the pendulum pivot bearings, the pendulum potentiometer, 
and the pendulum tachometer have been neglected. [Reference 
1] 


Differentiation of the above equations yields: 


m s(t) +mL (t) cos 9(t) - mL o*(t) sin 9(t) 

= in Ce) (6) 
See sin b(t) = ml o-(t) cos o(t) = V(t) -mg (7) 
J 3(t) = V(t) L sin 9(t) - H(t) L cos 3(t) (8) 
M s(t) = u(t) - H(t) - F s(t) (9) 


Note: For simplification, all independent variables which are 


functions of time will have (t) assumed but not displayed. 


ES 





Equation Limé@arizaetion 


Assume that the mass of the pendulum (m) is small with 
respect to the mass of the cart (M) so that H(t) may be ne- 
glected on the cart motion. 

Assume that the perturbation angle of the pendulum will 
always be small so that the small angle approximation may be 
used to linearize the equations. Therefore, sin ¢9(t) ~ o(t) 
ema cos 9 (t) - 1. Also neglect products of higher order 
derivatives of 9. 

By combining the pendulum and cart equations and then 


rearranging them, one obtains: 


— ae i 
so= = 8 + (10) 
-mL L 
6 = —S5s+ —— 39 (11) 
gms J +mL 
Tee ; ’ 
Let L' = a and consider it to be an "effective pendu- 
lum length." The equations of motion are now: 
_ oF : i 
Ss o= ast yu ao 
S _ Peele ace 
0 ~ mae a 2 (13) 


The equations of motion have been significantly simplified 
through the linearization process. It should be noted how- 
ever, that the actual physical inverted pendulum must oper- 


ate within a narrow limit of pendulum angle, be constructed 
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with very little real friction, and use a light weight, 
slender pendulum rod to actually behave as the mathematical 


model predicts. 


Seeley eat LON OF SYSTEM PHYSICAL CONSTANTS 

In designing the cart, values needed to be determined for 
all the coefficients in equations of motion (12) and (13), for 
meaningful computer simulation to take place and for design 
of the actual controller. Initial computer work was conducted 
using "ballpark" values to get a feel for how the system would 
react when stabilized with various control laws applied. The 
primary unknown value which must be found is F, the viscous 
Meret iogn coefficient of the cart. 

From measurements of the actual cart and pendulum, the 


Following data was recorded: 


Gact 
M = 8.45 lb = 3.83 kg 
Pendulum 
mee= 12.5 02 = 0.354 kg 
woe eeoOo ins = 1.016 m 
o Z J+mb- 
= ae = ! = ae 
J = DD nee. , LE Lie, ie aL 
therefore, 
, 4 
L = zh = O67 7 1 


iD 


i 








By rearranging equation (12) one obtains: 
Ms = u-FS§ (14) 


In this form, the equation may be recognized as Newton's 


Second Law: 
ec y HOEEEe'S 


The terms on the right side of equation (14) represent 


the total force exerted on the cart. 


u- F : (ole 5i) 


Fh 
1th 


The force exerted on the cart by the motor and gear train 
through the drive wheels, u, is simply a linear function of 


the voltage applied to the motor. Therefore one can define 


(56) 


& 
TT 
th 
< 


where r. 1S an applied voltage force (proportionality) con- 


stant and V is the voltage applied to the motor. Now, 
f = £ V-Fs hele? ) 
V 


In order to determine values for t., and F, it was 
necessary to eliminate selectively one of the two terms from 
equation (17). 

With the cart restrained from moving, a voltage was applied 


to the motor until the wheels started to slip. The force 
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exerted by the cart was measured uSing a Catilian Force Gauge. 


Weeeh no motion, S$ = 0 and f = fv Or; 


(18) 


th 
T 
<a} 


A series of test runs was performed which yielded an 


average applied voltage force coefficient of 


Some G2 Mbr/volt = 1.2286 newton/volt 


th] 
< 


Rearranging equation (14) and including the new definition 


Greu results in 





> = a a (19) 
s s 


Values for average velocity Ss, and average acceleration 
S must now be determined using a specific voltage to drive the 
motor so that a value for F can be determined. 

The pendulum was replaced with an equivalent mass of steel. 
The cart was run on the floor to observe motor tachometer and 
wheel potentiometer output on a brush recorder as a function 
of time as the cart accelerated from rest. A constant voltage 
of 15 volts was applwedeto the motor. 

Using the wheel potentiometer recorded output (voltage vs. 
time), the potentiometer's measured characteristics (degrees 
vs. voltage) (from Appendix A), and the physical size of gears 


and wheel, the average velocity was determined to be 


> 


Ss = ¥V = = = 4128 meters/sec (20) 


C 


EF 





Usang @ samidver method, the output from the motor 
tachometer versus time was converted to an average accelera- 


teeny OL* 
ee 9.75 meberns/eec- ‘ (210) 
Now using equation (19), one may solve for F. 


F = 32.922 kg/sec 


The equations of motion with physical values included are: 


-9.589 s + .26l u (22) 


wn: 
li 


AiG, ee ee (ah 


© ; 
il 


These values are within a factor of 5 of the values ob- 
tained in Reference 2 and within a factor of 8 of the value 
used in Reference l. 

In order to more easily manipulate the equations for the 
inverted pendulum, it will be convenient to express them in 
state variable notation. Straight forward checks for con- 
trollability and observability are available using state 


variable methods. 


C. STATE VARIABLE FORMULATION 
Intuitively one recognizes that with no forces applied 
except gravity, the inverted pendulum is unstable. When the 


pendulum shaft is balanced in a vertical position with the 
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puvotewat the bottom of the Shaft, the slightest perturbation 
will cause the pendulum to tumble. 
In general, any system may be checked for stability by 
looking at the roots of the characteristic polynomial, 
det [SI-A] = 0, where the system equations are in the state 
variable form x = [A]x + [Saoeeuenea stable system, all 
characteristic roots will lie in the left-half (negative) com- 
plex plane. Any roots lying in the right-half complex plane 
are considered unstable poles and lead to an unstable system. 
The inverted pendulum system equations of motion may be 
put into state variable form by defining the following states 


[Reference 1]: 


x, = s(t) (24) 
Ky S (t) ay 
x, = s(t) + L'9(t) (2Z10)) 
oy = Nene + L9(t) aT 


where, 


He is displacement of the cart, 
ec 1s velocity of the cart, 


X. 1s a linearized approximation of the displace- 
ment of a point on the pendulum at a distance 
L' from the pivot. It can be considered as 
the displacement of the pendulum, 


Ba is the linearized approximation of the velocity 
of the pendulum at the same point as made up Xa. 


Me) 





Therefore, 


x, = X, (28) 
x5 = + Xx. + 3 u (29) 
x. =), (30) 
x, = = S 8 = ce (31) 


x = 0 1 0 Onliesc 0 u 632) 
-F i 
0 oar 0 0 m 
0 0 0 iv 0 
=e wen 
nm 0 mn 0 0 
Pwo sl oRMiNATION OF STABILITY 
The characteristic equation of x = {A]x + [B]u comes from 
det[SI-A] = 0, 
[SI-A] = S = 0 0 (a3) 
F 
0 S+a 0 ) 
0 0 5) -1 
J _ me 
il e igh : 
ie _ Seeteee 2 oo 2 FG) lk 
det{SI-A] = S(S aS “7s aD = (34) 





The roots of the characteristic equation give the eigenvalues 
of the system: 


= 2.0, Ss = — S. = =76/L', aa g/L! 


Using the physical values as determined for the actual system, 


the eigenvalues become: 


The eigenvalues S5 and S 3 are negative, therefore stable poles. 


Eigenvalue Si 1s considered neutrally stable, and S 4 1s an 
unstable pole. 
One of the objects of this thesis was to stabilize the 


unstable system using a suitable control law, so controlla- 


bility of the system must be determined. 


Econ LROLEABLILLITY 

To solve control problems, such as the inverted pendulun, 
it is important to know if the system has the property that 
it may be moved from any given initial state to any other given 
State in a finite period of time. Controllability implies 
the ability to move the poles of a system from some starting 
position to any arbitrary place in the left-half complex 
plane. 

A theorem on controllability (Reference 1) states that: 


The n-dimensional linear time-invariant system 


Zale 





is completely controllable if and only if the column vectors 


Gemene: controllability matrix 


Pics age)... ACB) (35) 





Span the n-dimensional space. 
A matrix spans the n-dimensional space and is therefore 


Gontrollabbe if 


det{P] # 0 


Using the system numerical values, the state matrix 


equation is: 


x = 0 1 0 Cine + O | uu (36) 
0 -8.589 0 0 261 
0 0 0 l 0 
| 24. 486 0 14.486 0 0 j 
Pp = 0 apie oe VOMPEMnl an 
0, Jone = 2.242 19.254 -165.374 
0 0 0 aeeil 


Q Q eons 32.474 | 


Leen 


det{P] = 0.974 


Therefore, the inverted pendulum system is controllable. 
Although a system may be completely controllable, there 


may be times when all states of the system may not be available 
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for measurement. In order to find a control law which 
stabilizes a particular system which has unmeasurable state 
variables, one is led to the concept of observability. 
Knowledge of where the system is must be available before the 


system can be controlled. 


te OBoOLRVABILITY 
The concept of observability 1s based on an ability to 


observe the output vector of a system. 


y = ([C]x + [Du (38) 


A system is said to be completely observable if every 
peeeercoteace X(0) Can be determined from the observation 
mee) Over a finite time interval. 

A theorem on Observability (Reference 1) states that the 


n-dimensional, linear time-invariant system 


[xs 
T 


eae et 
(39) 


ee ss 


is completely observable 1f and only 1f£ the row vectors of the 


observability matrix 


on C (40) 





} 9 
Pee 
| can" 


—! 
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span the n-dimensional space. The matrix spans the n- 
dimensional space if det[Q] # 0. 
With the fourth order system under consideration there 
are 15 possible output combinations for the vector [C]. 
When the output vector is made up of states xy and x 


eB 
meoeeckample, y = {111 0]x and: 


Q = 1 0 1 0 (41) 
0 1 0 1 
-14.486 -8.589 14.486 0 
0 59.285 0 14.486 | 
Beato = =1297.917 


Therefore the inverted pendulum is observable for the output 
y = [1010]x. 
Observability of the system for all output combinations 


of {C] is as follows: 


Observable Not Observabie 
Omen: (0 00 1] 

Oro) | LOMECr0) 

On Pero: Om InO «1 

Powel ise! | Tela) (8) 883 

Cleve 8 ie O50 } 

pres 0') 

(101 1] 


pee] 





Peo) 


paae 1) 


With the conditions of controllability and observability 
known for the system, the analysis and design of a linear 


control law for stabilizing the inverted pendulum may proceed. 
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IIIf. STATE VARIABLE FEEDBACK 


Pwo ltatet, VARLABLE FEEDBACK CONTROL INTRODUCTION 

A control system is a dynamic system which, over a period 
of time, behaves in a manner prescribed by the control law 
which describes the control system. In the case of an auto- 
matic control system, human input is not required and may 
not be desired. 

The primary sections of a control system are: the plant, 
which is the system to be controlled; one or more sensors to 
provide information about the system; and a controller which 
compares measured values to desired values and changes plant 
inputs accordingly to arrive at the proper output. 

One of the functions of the controller is to move the poles 
of the plant to locations in the left-half complex plane so 
that performance of the system may be improved. When the 
plant by itself is unstable, the main function of the con- 
troller is to stabilize the system by moving the closed-loop 
poles to proper locations in the left-half plane. 

A feedback control system is one which attempts to main- 
tain a set relationship between the output and some reference 
input by comparing them and using the difference signal as a 
means of control. 

State variable feedback is particularly convenient because 
ew oeeee<ecONtains all pertinent information about the sys- 


tem. The basic restraint placed on the analysis in this section 
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is that the complete state of the plant can be measured 
accurately at all times and is available for feedback. 
For a linear time-invariant system such as x = [A]x 


ei, 2 time-invariant control law which applies is: 
u = -({K]x + u' (42) 


where [K] is the feedback matrix and u' is some input such 
as from a human operator. In the inverted pendulum system 
ees anored. 


The control law 
u = -K,x, - K.x. - K Xq = i 2x (43) 


applies for a fourth order system. 
Substitution of this control law into the system state 


niet (32) yields: 


= 0 1 0 0 x (44) 

ee aE eke. digi 

M M M M M 

0 0 0 1 

air E'S 

at : im : 
Oe 

= ene Epc = ee er) 
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ug 
tes 
Kee K 
4 2 5 ees Z Boa _ _ 
(46) 
Pom a = 
~ rin | Ky K) 0 


With system numerical values included, the equation becomes: 


x 


Pee (ye) K.+ 8.589)S> + (.261 K. - 14.486)s° 


Z 1 


(47) 


+ 3.782(-K, ee ater) ct 3.782 (-K, aad 


2 3) 


With proper selection of feedback gains one may arbitrarily 
assign pole locations to stabilize a system or to simply 


improve its performance. 


B. FULL STATE VARIABLE FEEDBACK CONTROL 

Another method for determining feedback coefficients for 
a time-invariant control law involves the use of "engineering 
judgment" to place the poles of the closed-loop system in the 
left-half complex plane at a position of one's choice. By 
choosing the poles far to the left in the complex plane, the 
convergence to the zero state can be made arbitrarily fast. 


To make the system move fast however, requires large input 





amplitudes, and in an actual physical system, there is a 
limit as to how large the input can be. 


Asisume it is desired that the poles be placed at 


The desired characteristic equation is then 


PRECaeNG. + 6235S tulle? Sitene2 = oO». (48) 
The control law which applies in this case is the same 


as for the general state variable feedback case 


Again for the system to be stable, each coefficient of 
the characteristic equation must be positive. Equating 
coefficients of the augmented characteristic equation (47) 
and the desired characteristic equation (48) yields the proper 


feedback coefficients for pole placement 


Ky = 296.9 

Ko = on wo 
K Se ie 
Ky = -88.825 


Computer simulation using the Interactive Ordinary Differ- 
ential Equation Solver (IODE) on the Naval Postgraduate School 


IBM 370 shows that pendulum stabilization to a displacement of 


ao 


1) 





.007 meters in 2 seconds should take place from an initial 
condition of .0667 meters (.1 radians) of pendulum angle. 
(See Appendix B for simulation results.) 

Figure 2 is a block diagram representing the unstable 
plant of the inverted pendulum and the state variable feed- 
back controller which stabilizes it. 

Full state variable feedback has been shown to work well 
for a system with all states accurately measured. 

The actual cart does not respond in the same manner as the 
Simulations would predict because the cart 1S receiving a 
continuous input from all of the sensors instead of a single 
initial disturbance. As might be expected, once the actual 
pendulum system motion is started, the pendulum angle oscillates 


on either side of zero. 
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ye eee tet ekUClULON OF INVERTED PENDULUM MODEL 


Building a working model of the inverted pendulum for this 
thesis was quite a challenge. Little technical data could be 
obtained concerning construction of any previous cart and 
pendulum systems, therefore design ideas came from line draw- 
ings, photographs of other systems, and basic intuition. 

As a starting point, the system was to be portable for 
viewing in a classroom environment rather than as a laboratory 
setup. Large equipment such as analog computer and heavy 
power supply could not be used, therefore a solid state power 
supply and controller were built. Target weight for the 
cart assembly was arbitrarily set at 6 pounds, and a pendulum 
length of 40 inches was chosen. 

The motor/generator, pendulum tachometer, pendulum potenti- 
ometer, and wheel potentiometer were chosen primarily because 
they were already available within the N.P.S. Aeronautical 
Engineering Department. The physical dimensions of the cart 
came about simply from the physical sizes of the motor and 
sensor components to be used. Rubber wheels were chosen to 
give good traction and to minimize damage to table or bench tops 
when the cart was demonstrated. A power supply and control 
Circuit (covered later) were housed in a small carrying case 
and an umbilical cord was used to supply power and control 
Signals to the cart and receive sensor feedback from the cart 


and pendulum. 
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Photographs of the cart assembly and control box are in 
pueures 3 through 6. 

Testing and measurement of the tachometers and potenti- 
ometers was carried out in the Aeronautical Engineering Elec- 
tronics Laboratory. Calculations, tables, and plots of the 
sensor constants are contained in Appendix A. 

A full state variable feedback control law was developed 


in Section III which placed the closed-loop poles at 


The control law was 


u= “Ki xy = K5X.5 = KX - KyX4 (49) 
where 
Ky = 296.90 
Ko = diesen: 7 
K 4 =) =-359 72 
Ky = -86. de 
Recall that 
x, = 5 oe = Cart vosition 
Xx, = S S = cart Feleocliy 
wa = ~=6O St Ey! > = pendulum position 
x, = S + ty b = pendulum velocity 
ie = 0.677 m L' = effective pendulum length. 
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CART QUARTER VIEW 


Figure 3. 
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Substitute state definitions into eguation (49) to get 


u = (-K, -K,)S f: (-K,-K,)S = K,L"9 aa K,L'¢ : (50) 


Equation (50) is the control law as referenced to the 
measurement variables of the system constructed. With the 


gain constants substituted, equation (50) is 
ee ees 76S + 230.0 6 + 60.14 4 . (51) 


In order to implement the control law, each of the measured 
variables must be referenced to a common base, amplified 
appropriately, and summed with the other variables to produce 
ame centrol. 

From Appendix A, each sensor was measured, and a constant 
was developed which references each sensor to a base level of 
l1 volt. In the measurements, the range of operation of each 
sensor was taken into account in weighting the constant. 


The final values obtained for the sensors used were: 


K = .1734 meters/volt --- wheel potentiometer (S) 

K = .2012 meters/second-volt. --- motor tachometer (S) 
= = 1.1073 radians/volt --- pendulum potentiometer (9) 
Kn = 9.39 radians/second-volt. --- pendulum tachometer (4) 


An additional potentiometer was included in each sensor 


input circuit which acts as a voltage divider. The voltage 
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divider actions provide signals from each sensor which are 
of the proper magnitude before being amplified and summed 
for the control law. 

The sensor calibration constants must be applied to the 


control law to realize a properly proportioned electronic 


G@ntroller. 
uo= K (7K, -K)s + K (7%, 7K,) Ss - at o- KU K,L d (52) 
or for the system desired, 
u = K (42.85)S > K 6 79-76)8 + ae $ 
+ K_ (60.14) » (3) 


and finally, after substitution of actual sensor constants 
u = 7.43 § + 15.24 S + 254.68 9» + 564.66 4 (54) 


Implementation of the controller using operational ampli- 
fiers and power transistors is shown in figure 7. 

Each front panel potentiometer has a range scale from 0 
to 10 using a calibrated dial. 

The gain adjustment potentiometers size the measured 


variables as follows: 


Potentiometer one S 
Potentiometer two S 
Potentiometer three b 


Potentiometer four $ 
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G@ieweernal Controller circuitry provides an additional 
amplification to supplement the voltage division of the gain 
adjustment potentiometers. An infinite number of pole place- 
ment positions are available through variations in the state 
and can be 


variable feedback gain constants K, through K 


1 4’ 
implemented with the constructed controller. The only limi- 
tation is on the maximum coefficient values in equation (54). 


The values which may be used are: 


0 to 100 S 
0 to 100 S 
OO 0100 b 
0 to 1000 4 


aewsece up the controller for the control law of equation 
(54) dial in the following gain adjustment potentiometer 


Seeerngs on the control box. 


pot #1 .74 (id SS) 
pot #2 ies (15.24 $) 
pot #3 2.55 (254.68 9) 
pot #4 5.65 (564.66 9) 


Potentiometer number 5 adjusts the overall magnitude of 
PienOuEpUEsVOlEAGe £O Ehe cart drive motor. It was found that 
a value of 5.0 on that potentiometer worked well. Potenti- 
ometer number 6 is not used at this time but was installed 


for future expansion. 
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By experiment, the cart works best with the following 


settings: 
pot #1 SD 
pot #2 O (or a small value) 
pot #3 10 
pot #4 10 
pot #5 SU 


The inverted pendulum system as constructed does not work 
as well as hoped. Ability to closely control the cart's 
position has thus far been unsuccessful. The cart's region 
of motion spans approximately 8 to 12 inches of horizontal 
cart travel. Cart position oscillates on either side of a 
zero location with an increasing divergence and an increase 
in the signal from the wheel potentiometer only causes the 
system to go unstable more rapidly. When values above zero 
are input from the motor tachometer, the cart motion becomes 
more sluggish. Best cart performance and longest stability 
1s achieved with the motor tachometer output grounded (potenti- 
ometer set to zero). The controller noticeably attempts to 
correct the cart position after it has travelled 4 to 6 inches 
from center, however, the motor inputs are somewhat large, 
and after the cart moves 8 to 10 inches from center, the motor 
response appears to be too large and the pendulum swings out 
of the stable range and topples. At this time there is not a 


full understanding of the problem with cart velocity. 
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The range of operable pendulum angle is small. In order 
to have the system work, one must help the pendulum achieve 
a small stable oscillation before letting the system to on 
its own. The system operates for 10 to 15 seconds then goes 
unstable with larger and larger movements in both cart and 
pendulum. This observation can only lead to the conclusion 
that the system is only neutrally or marginally stable at best. 
Because the motion of the cart seems to be driving the pendu- 
lum unstable, it is possible that the pendulum poles are ade- 
quately positioned on the left-half complex plane and the cart 
poles are on the axis or even in the right-half plane. 

There is a sizeable amount of real friction in the pendulum 
mount from bearings, potentiometer and tachometer drag which 
have been ignored in the development of the controller, but 
which possibly contribute to the instability of the system. 

Some recommendations for areas of possible improvement 


are covered in the final section of this paper. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


i  CGONCLUSIONS 

The purpose of this paper was to design an inverted pendu- 
lum which could be used to compare theory with experimental 
results. In that sense, the purpose was not met. 

Equations of motion were developed for the inverted pendu- 
lum and converted into state variable form. State variable 
feedback was selected as a basis for control of the actual 
system constructed. In order to reduce the complexity of the 
system for control analysis, the equations of motion were 
linearized. Friction in the pendulum pivot was ignored, but 
in the actual system, friction was definitely present. 

Friction possibly added damping to the system so that the 
pendulum fell more slowly, however, it would have also forced 
the cart to make larger and longer duration corrections to 
Zhen ure pendulum to a vertical position. 

The size of the cart drive motor along with the traction 
of the drive wheels could possibly be the cause of the system's 
inability to become stable. 

It was assumed that the dynamic response of the motor 
and its electronic Circuit was sufficiently fast to be con- 
Sidered instantaneous; that switching from forward to reverse 
rotation was done without delay. 

It is speculated that the actual motor and gear drive 


reached their acceleration limits when trying to catch a 
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pendulum angle much above +5 degrees. The motor simply could 


not accelerate the cart fast enough to catch the pendulum. 


B. RECOMMENDATIONS 

The following list includes several possible ways to 
improve the physical system. 

1. Reduce friction in the pendulum pivot by: using a 
small, low friction, high precision potentiometer or magnetic 
pickup device to sense pendulum angle; replacing the permanent 
magnet pendulum tachometer with a magnetic drag-cup tachometer 
or with a micro-lightweight accelerometer attached to the 
pendulum at its center-of-gravity or tip. 

2. Purchase a high quality motor/tachometer assembly for 
driving the cart. Choose a motor with a torque output in the 
50 to 100 ounce-inch range and be certain that all constants 
are available for the new motor and tachometer so that a full 
analysis may be made using the appropriate constants. 

3. Install tires which are at least twice as wide and 
perhaps twice as tall as those now used to provide more con- 
tact area for improved traction. 

4. Analyze the system using sensor and motor constants 
(if available), gear ratios, moments of inertia of all rotating 
pacts and friction. 

5. Study the effects of measurement errors, system dis- 
turbances and noise on observers along with optimal methods 


for selection of gain matrices for pole placement. 
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6. Implement cart controllers using full-order and 
reduced-order observers to drive the controller. 

7. Study the effects on stability of using slightly 
flexible pendulums and pendulums of various lengths. 

8. Study the possibility of implementing the control 


system digitally using a microprocessor chip. 
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APPENDIX A 


SENSOR MEASUREMENTS 


Measurements for determination of sensor constants was 
accomplished using equipment in the Naval Postgraduate School 
Aeronautical Engineering Department Electronics shop. 

Al. Pendulum potentiometer measurement 

A2. Wheel potentiometer measurement 

A3. Motor tachometer measurement 


A4. Pendulum tachometer measurement 
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Al. PENDULUM POTENTIOMETER MEASUREMENT 


RAW ADJUSTED VOLTAGE 
REFERENCE OUT UL 
(DEGREES) (VOLTS) 


REFERENCE 
(DEGREES ) 


28 

30 

50 

70 
90 
oo 
130 
> 0 
ie 
9'0 
210 
230 
20 
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290 
5) 1E8 
S610 
Bi 
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2 
22 
42 
62 
82 
OZ 
eZ 2 
142 
or 
SZ 
ZU2 
ee 
242 
Zoe 
262 
B02 
Bi 
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TABLE 1. PENDULUM POTENTIOMETER MEASUREMENT 


From the slope of the measurement curve 


volts a ony 
degrees  285.5° 


Ale 


1 volt corresponds to 
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.0315 V/degree = 


Poo Vokes, radian 





eoO SAL 


= 1.1073 radians = Kp 





For operation on the actual pendulum, the potentiometer 
was operated at the center of its range with +15V applied to 
one end and -15V applied to the other end, giving a zero 
volt reference in the center at 163.5°. Operational range 
of the potentiometer is indicated on the measurement curve 


figure 8 and reflects a pendulum angle of +8 degrees. 


A2. WHEEL POTENTIOMETER MEASUREMENT 








pot voltage voltage 
earn output output 
number (volts) (volts) 
a 
l 
0 os 
90 1.88 
“ 180 2.25 7 
270 Dee 
0 3.01 0 10.50 
: 90 3.37 : 90 One 7 
180 a) 9G | 180 i 2 5 
a7 0 13 270 LiL ae 
0 Aeash ( 00 
; 90 4.89 i 90 12138 
180 5.26 180 76 
270 ace | 270 pee 2 
| 0 amon 0 | 13.50 
: 90 638 - 90 (13.88 
| 180 575 180 )14.25 
a 70 Te 270 114.63 
360 115.00 


— 


TABLE 2. WHEEL POTENTIOMETER MEASUREMENT 
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From the slope of the measurement curve, figure 9: 


Oo o>. OL volts 


720 degrees Se aes 


igfemcarte travels 38.7 cm for each revolution of the 10- 


turn wheel potentiometer, therefore: 


“- Beer wcu, tamecer , 1 deg 


: mecca oomem sages volt ~ -+/3% ™/volt 


1 volt corresponds to .1734 meter. 


A3. MOTOR TACHOMETER MEASUREMENT 





















POWER 
SueP LY 
OUTPUT 
(VOLTS) 


TACHOMETER 
OUVe UT 
VOLTAGE 

(VOLTS) 


JAGABI 
INDICATOR 
(RPM) 


STROBE 
LIGHT 
(RPM) 


386 

S15 
1340 
350 
Ze 0 
20.90 
3420 


3920 


iheor 3. “MOTOR TACHOMETER MEASUREMENTS 
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From slope of the measurement curve 


18 -2 V a om eo ve 023 V-sec 
3720 -430 RPM ~ 3320 RPM 100 RPM : rad. 


The cart travels 5.9 cm for each revolution of the motor/ 


tachometer. Therefore 


ge = o:9 CM x cae beWeRey , 1 min 
W Rev 130 sem ~4819 V-min 60 Sec 
= meters 
ae volt-sec 


1 volt corresponds to .2012 meters/second. 


A4. PENDULUM TACHOMETER MEASUREMENT 


Voltage 
Output 
(Dee. volts) 





TABLE 4. PENDULUM TACHCMETER MEASUREMENT 


Only three values were obtained for the pendulum tachometer 
as the only method available for turning the tachometer was 
uSing a three speed drill press then measuring RPM with a 


strobe light. 
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eee, _ 


- _100 Rev_ , 1 Min , 27 Rad , 360° _ 4g 39 _radians 
An 2.21V -min 60 Sec 180° Rev E volt-sec 


1 volt corresponds to 9.39 radians/second 
Data plate information on this tachometer specifies a 


2.1 volt/100 RPM linear response which is very close to the 


experimentally derived value of 2.21 V/100 RPM. 
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APPENDIX B 


PROBLEM SIMULATION 


Computer simulation was done using the Interactive 
Ordinary Differential Equation Solver Program (IODE) on the 
Naval Postgraduate School IBM 370 system. 

All simulations use actual physical parameters for the 
inverted pendulum system which was constructed. Simulations 
include tabular output as well as french curve fit plots 


from the Hewlett Packard HP-9830 graphics system. 
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